Gap nodes induced by coexistence with antiferromagnetism in iron-based 

superconductors 
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We investigate the pairing in iron pnictides in the coexistence phase, which displays both su- 
perconducting and antiferromagnetic orders. By solving the pairing problem on the Fermi surface 
reconstructed by long-range magnetic order, we find that the pairing interaction necessarily becomes 
angle-dependent, even if it was isotropic in the paramagnetic phase, which results in an angular vari- 
ation of the superconducting gap along the Fermi surfaces. We find that the gap has no nodes for 
a small antiferromagnetic order parameter M, but may develop accidental nodes for intermediate 
values of M, when one pair of the reconstructed Fermi surface pockets disappear. For even larger 
M, when the other pair of reconstructed Fermi pockets is gapped by long-range magnetic order, 
superconductivity still exists, but the quasiparticle spectrum becomes nodeless again. We also show 
that the application of an external magnetic field facilitates the formation of nodes. We argue 
that this mechanism for a nodeless-nodal-nodeless transition explains recent thermal conductivity 
measurements of hole-doped (Bai_ a! K a ;)Fe2As2 [J-Ph. Read et al arXiv:1105.2232 . 

PACS numbers: 



I. INTRODUCTION 

Exploring the different regions in the phase diagram of 
iron based superconductors (FeSCs) is an important step 
towards developing a unified understanding of the physics 
of superconductivity in these systems. A typical phase 
diagram of a FeSC in the (T, x) plane, where x is dop- 
ing, shows a metallic antiferromagnetic order, also called 
spin density wave (SDW), below TV at X = 0. Upon 
doping, the SDW order parameter is suppressed and su- 
perconductivity (SC) emerges with maximum T c (x) near 
the point where T c (x) exceeds Tjv(x) (an "optimal dop- 
ing"). The gap symmetry at optimal doping is most 
likely s^ , but the structure of the s-wave gap varies 
from one material to another: the gap is nodeless in 
Ba(Fei_ x Co x ) 2 As^- and (Bai_ x K a; )Fe2As^-, but has 
nodes in BaFe2(Asi_ x P a; )2^ 1 ^. 

In this article we investigate the behavior of 
the s^ SC gap in the underdoped regime of 
the Ba(Fei_ 2; Co a ;)2As2 and (Bai_ x Ka;)Fe2As2 systems, 
where T c (x) < T/v(x), and superconductivity emerges 
in a continuous fashion from a pre-existing SDW or- 
der. The microscopic coexistence of SC and SDW be- 
low T c (x) has been reported by magnetization and NMR 




nodeless < 



nodeless 



experiments^ — . 

The electronic structure of FeSCs in the paramag- 
netic phase consists of near-circular hole pockets and el- 
liptical electron pockets. We assume, following earlier 
worksi 3 - - — , that the interaction leading to SDW order is 
independent on the angle along the FS 20 . In this situa- 
tion, the Brillouin Zone (BZ) is reduced once the SDW 
order sets in, and the Fermi surface (FS) is reconstructed 
into four banana-like pockets (right panel in Fig. []}. As 
the SDW order parameter M increases, the reconstructed 
pockets shrink and eventually disappear. This happens 



Figure 1: Schematic phase diagram showing the nodeless- 
nodal-nodeless transitions inside the region where supercon- 
ductivity (SC) and spin density wave (SDW) coexist. Each 
transition is roughly associated with the disappearance of one 
pair of magnetically reconstructed Fermi surface pockets. Mi 
and M% refer to the values of the SDW order parameter M at 
which the two pairs of pockets disappear respectively. A dip 
in T c is observed near M2, when both pockets disappear. 



in two stages: first, at Al = Mi, the number of pock- 
ets shrinks from four to two, (middle panel in Fig. [I}, 
and then the remaining two pockets vanish at a larger 
value M — M% (left panel in Fig. [TJ. Because these 
reconstructed FSs are formed by mixing electron and 
hole bands, on which the s" 1 SC order parameter has 
different signs in the absence of SDW order, one could 



naively expect that the SC gap on the reconstructed FSs 
alternates between plus and minus signs, implying that it 
must have nodes. However, calculations show 13 that this 
is not the case because SDW order, with e.g. (S x ) 7^ 0, 
mixes hole and electron FSs with opposite a z compo- 
nents of electron spins. Because a spin-singlet s-wave 
gap changes sign under a z — » —a z , the sign change im- 
posed by the s H gap structure is compensated by the 
sign change due to the flip of er 2 , and, as a result, the 
gaps on the reconstructed FSs retain the same nodeless 
structure as in the absence of magnetic order. 

Experiments near the onset of the co-existence phase, 
where M is small, are in agreement with this reasoning. 
For instance, the ratio k/T, where k is the in-plane ther- 
mal conductivity, tends to zero in the T — > limit, as 
it is expected for a superconductor with a nodeless gap. 
However, recent thermal conductivity measurements on 
(Bai_ a; K a ;)Fe2As2 deep in the coexistence phase 22 found 
that there is a range of dopings where k/T remains finite 
at T = 0, as it happens in a superconductor with line 
nodes. At even smaller dopings, k/T again vanishes at 
T = 0. Application of a magnetic field makes the doping 
dependence more smooth and extends the range of nodal 
behavior. 

In this article we explain the sequence of nodeless- 
nodal- nodeless behavior as a transition from a node- 
less s^ gap in the region where all four reconstructed 
FS pockets are present, to an s~* gap with accidental 
nodes in the region where two FS pockets are present, 
to the fully gapped quasiparticle spectrum in the region 
where the remaining two reconstructed FS pockets are 
also gapped (see Fig. [1}. We show that the pairing in- 
teraction in the coexistence phase acquires an additional 
angular dependence as it gets dressed by angle-dependent 
coherence factors associated with the SDW order. The 
angular dependence of the interactions in turn gives rise 
to an angular dependence of the superconducting gap. 
This effect is weak in the limit of small M , considered in 
Ref. Il3l . when all four reconstructed Fermi pockets are 
present, but is strong and gives rise to gap nodes when M 
becomes large enough to gap out one pair of Fermi pock- 
ets. At even larger M, deeper in the coexistence region, 
the remaining pair of Fermi pockets is gapped out. In 
this situation, superconductivity still develops over some 
range of dopings 14 , but the quasiparticle excitation spec- 
trum remains gapped irrespective of the structure of the 
SC gap, and k/T again vanishes at T = 0. Concurrently, 
the superconducting transition temperature T c shows a 
minimum when the remaining pockets vanish. The the- 
oretical behavior across the SDW-SC coexistence region 
is shown schematically in Fig. [1] It agrees well with the 
experimental results of Refs. l22l ? 

The paper is organized as follows. In Sec. HI we de- 
scribe the formalism listing out the details of the FS ge- 
ometries, the nature of the reconstructed FSs, and the 
generic form of the SC gap in the coexistence region. In 
Sec. IIIII we discuss the solutions for the gap obtained 
within this formalism, show when and how nodes arise 



in the coexistence region, and provide analytical expla- 
nation of the results. We also briefly discuss the role 
played by bands that do not participate in the formation 
of SDW. In Sec lIVI we show that the doping range where 
the SC gap has nodes is enhanced by applying a magnetic 
field. We state our final conclusions in Sec. \V\ 



II. THE PAIRING PROBLEM ON THE 
RECONSTRUCTED FERMI SURFACE 

Our point of departure is a microscopic band model of 
interacting fermions located near hole or electron pock- 
ets. The electronic structure of FeSCs consists of two or 
three hole pockets, centered at (0, 0) in the folded zone 
(two Fe atoms per unit cell) , and two hybridized electron 
pockets centered at (71", n). The interactions between low- 
energy fermions are generally angle-dependent already in 
the normal state, before either SDW or SC order sets in, 
which gives rise to angle-dependent SDW and SC gaps 
even outside the coexistence region. 

As we said, we follow earlier works and consider a sim- 
plified model in which we neglect the angular dependen- 
cies of the interactions in the normal state, and approxi- 
mate the normal-state interactions as constants^ 3 . - — . We 
further neglect the hybridization and perform calcula- 
tions in the unfolded zone (one Fe atom per unit cell) , in 
which electron pockets are ellipses centered at (0, 7r) and 
(71-, 0) and hole pockets are at (0, 0) and (it, tt). Finally, we 
neglect the third hole pocket and only consider two cir- 
cular hole pockets centered at (0, 0) in the unfolded zone 
(in orbital notations these are d xz /d yz pockets). Earlier 
works have foundr-^ that SDW magnetic order emerges in 
this model with ordering vector Q = (0, ir) or Q = (n, 0), 
mixing one hole and one electron pocket. For definiteness 
we set Q = (0, n) in which case one hole pocket at (0, 0) 
and one electron pocket at (0, n) are mixed up. The other 
two pockets are not participating in SDW and remain in- 
tact once the SDW order sets in. In line with our goal, 
we first consider only the 2-pocket model with one cir- 
cular hole and one elliptical electron pocket, which get 
reconstructed in the SDW phase. We discuss the role of 
the other pockets later in the next section. 

To quadratic order in the fermions, the Hamiltonian 
of the 2-pocket model in the SDW phase is Hq + Hsdw, 
where 



Hn = 



7 , (^k c ker C k<x + e k+Q/k+Qcr/k+Q<T 
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DW 
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Here c and / are fermionic operators for hole and elec- 
tron states, respectively, a is ±1, and M is the SDW 
order parameter, which we treat below as a variable. 
In practice, larger M correspond to smaller dopings, 
whereas smaller M correspond to doping close to the op- 
timal one. The fermionic dispersions are 
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where to, m x , rn y are the effective band masses of the 
fermions. For definiteness, we set the interatomic spacing 
to one and use jjl c = \hf = /x, m = l/(2/i), m x — 0.5/(2/i), 
and TOj, = 1.5/(2/x). These parameters are chosen to give 
the FS geometry as in Fig. [5] We will also measure M 
in units of 2fi. 

The M term in fTJ) couples c and / operators such that 
the eigenstates of (fl} are coherent superpositions of elec- 
trons and holes. These are described by new fermionic 
operators a and 6, with dispersion £"£' that vanishes at 
the reconstructed FSs. The transformation to the new 
operators is 



are the quasi-particle excitation energies of the SDW 
state (the plus sign corresponds to a fermions, the mi- 
nus sign corresponds to b fermions). 

The condition E^' = defines the new reconstructed 
FSs. For small enough M , there are two pairs of banana- 
shaped reconstructed pockets, the a-pockets and the b- 
pockets (see Fig. [2J, with the latter larger than the for- 
mer. At the critical M = Mi given by 



Mi = 



1 |to 



4 J 



(9) 



the a-pockets disappear while the 6-pockets remain (Fig. 
fright). At even larger M — M 2 given by 
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1 |m — to t 
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where 
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such that 



cos 20 k = 



sin 20k = — - 



M 2 + ( e -^' 



M 
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In terms of the new operators, 

Hq + HsDW = 2_^ ( ^k a k(T a ktT + -^k^kcr^kcr 

where 
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k.cr 
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+ M 2 (8) 



the 6-pockets also disappear, i.e. all electronic states 
become gapped by SDW—. For our set of parameters 
Mi k, 0.102 and M 2 « 0.177. 

We now turn to the issue of the pairing on the re- 
constructed FS. There are four residual interactions be- 
tween the original fermions located near the hole and 
electron FSs — : the interaction between fermionic den- 
sities near hole and electron pockets C/i cJ.c CT f*, f a i , 
the exchange interaction between hole and electron 
pockets C^2cJ./ (T /J,c (T ' , the umklapp pair-hopping in- 
teraction in which two fermions near a hole pocket 
are converted into two fermions near an electron 
pocket and vice- versa C/ 3 /2(4/ ff 4,/ ff ' + /Jc -/ ( t'<v) , 
and the density-density interaction within each pocket 
174/2(4 c CT c+,cv + ftUflfa')] ■ In the absence of SDW 
order, only U3 and U4, contribute to the pairing channel 
(C/3 must be larger than U4 for s H pairing). In the 
SDW phase, c and / operators are mixed up, and all 
four interactions contribute to pairing vertices for a and 
b fermions. The interactions U2 and U4 do not give rise 
to new physics and only renormalize the values of the 
pairing vertices obtained from the U\ and U3 terms. To 
shorten the formulas, we neglect the E/ 2 an d U4 terms 
and approximate Hi nt by 



Hint = ^2 ^l4 ./J <7 y/3tr'C4 

[1234] 



~Y \ C L4a'h<r J fa<r+ h -' 



(11) 



where X) [12341 denotes the sum over all momenta subject 
to ki + k 2 = k3 + k4, and the summation over repeated 
spin indices is implied. Converting to new operators via 
Eq. ([3]) and keeping only the interactions which con- 
tribute to the pairing, we obtain from Eq. pip 
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Figure 2: Reconstructed a and 6 Fermi pockets for several values of the SDW order parameter M. The a-pockets are smaller and 
disappear first upon increasing M. At even larger M, the 6-pockets also disappear (not shown). Each pocket was discretized 
into a finite number of points indicated in the left figure for one of the b pockets. Our convention is such that the numbering 
for each pocket starts at the tip of a "banana" and goes first along the "inner" side of the pocket and then along the "outer" 
side. 
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Note, however, that the three angular components of 
A a ^ are not orthogonal, since e.g. J dq\\ cos 20 q does 
not vanish. Substituting these expressions into the gap 
equation, we re-express it as a matrix equation for T c and 
9i- 



(12) 



where 

jjaa 



u. 



q ,k - u qk = -f (-1 + cos 26>q cos 26> k + s sin 2# q sin 26> k 
C/q 6 k = C/q a k = - -j- (1 + cos 26> q cos 26> k + s sin 26» q sin 26> k ) 

and we denned s = U1/U3. 

We see from fp~2|) that both intra-pocket and inter- 
pocket pairing interactions acquire angular dependence 
via the cos2# q , cos 26^ and sin 2# q , sin2# k factors. Con- 
structing the set of linearized BCS gaps equations for 
the angle-dependent gaps A a (q) and A&(q) by standard 
means we find at T = T r 
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where L = log ( ■£- ) , A is the upper cutoff of the theory, 

and the Fermi velocity vf = Uf^i) varies along a— and 
b Fermi surfaces. 

A generic solution of this gap equation is of the form 



A Q (q) = 91 + 92 cos 26» q + 53 sin 26» q 
Ab(q) = ,9i - gi cos 26» q - g 3 sin 26» q 



(13) 



Without SDW order, sin 29k — 0, cos 29k = sign(e c — 
e/), and N s , N sc , and N ss all vanish. Then g 3 = 0, 
while g 2 / gi = -N C /(N + ^ N 2 - N 2 ). As a result, 
the gaps on the inner and outer sides of the a and 
b pockets reduce to g x (1 + N C /(N + ^/N 2 - N 2 )) and 

g 1 (l-N c /(N+ y/N 2 - TV 2 )). These two coincide with 
the gaps on the c and /—Fermi surfaces at T c in the para- 
magnetic phase, which we set in our model to be angle- 
independent (we recall that the sign of the two gaps is 
the same because we attributed an extra minus sign to 
the gap on the /—Fermi surface by nipping the spins 
of /— fermions). The magnitudes of the gaps on c and 
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Figure 3: The quantities Nj/N, where Nj = N C ,N S ,N CC , 
N ss , plotted as functions of M. The dashed line indicates the 
value M — Mi at which the a-pockets disappear. The ellipsis 
indicate the area around Mi where the ratios rapidly evolve. 
Note that N c /N remains small for all M, while N s /N and 
N as increase with M and become positive and of order one at 
M > Mi. 



of the bananas. Then N s , N sc , and N ss are small, and 
N cc /N rj 1. Solving Eq. (0 in this limit, we find that 
.93 "C 32 <C g\ , hence A and A„ remain almost constant 
along the a— and b— pockets, including the tips of the 
bananas. 
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Figure 4: The behavior of cos 2#k (left column) and sin 2#k 
(right column) along the a and b pockets for M = 0.04 and 
M = 0.09. The gap is given by Eq. (13]). The "length along 
the pocket" refers to the numbers along each pocket specified 
in Fig. [3 



/— FSs are generally unequal as N c ^ because of the 
different geometry of hole and electron pockets. Numer- 
ically, however, for our set of parameters N c /N <C 1, 
i.e., at M = the gaps on a and b FSs are essentially 
isotropic. Once SDW order sets in, all Ni and Nij become 
non-zero, 53 becomes finite, cos 20k and sh^^k become 
smooth functions along the pockets, and the gaps A a (q) 
and Af,(q) acquire smooth angular variations. The mag- 
nitudes of the variations depend on the two parameters: 
the ratios 32/51 and 33/31 and the actual variation of 
cos2#k and sin26*k along a— and 6— Fermi surfaces. 

In the small M limit, considered previously in Ref. U3 , 
sin 2#k is small except for narrow ranges near the tips 




Figure 5: Top panel: the ratios (/2/31 (blue line) and gz/gi 
(red line) as function of s for M — 0.09 (solid lines) and 
M — 0.12 (dashed lines). These values are slightly below and 
slightly above the critical Mi = 0.102 at which two out of four 
pockets disappear. Both ratios flip sign between M = 0.09 
and M = 0.12. The key result here is that for s > 1, the 
magnitude of 33/gi becomes larger than 1 for M = 0.12. This 
gives rise to accidental nodes of the gap on the remaining FS 
pocket (see text). Bottom panel: the plot of the same ratios 
as functions of M for s — 5. The magnitude of 32/gi remains 
relatively small, while the magnitude of 33/31 becomes greater 
than 1 for M > Mi. 



For arbitrary M, 32/31 and 33/31 are rather complex 
functions of the ratios Nj/N and also s. We plot Nj/N in 
Fig [31 We see that N c /N and N cs /N, are small because 
the integrands for N c and N cs contain cos20fc, which 
changes sign between the inner and outer sides of the 
a and 6-pockets (Fig. 3]). At the same time, the mag- 
nitudes of N s /N and N ss /N increase with increasing M 
and become of order one. 

In Fig. [U we plot cos 2#k and sin 26>k along the a and 



6-pockets. We recall that the prefactors of the cos 26>k 
terms in A a and Af, are already small because 32/31 is 
small. Fig. [4] shows that the variation of cos 2#k along 
each of the pockets decreases with increasing M, making 
the g-2 cos 2#k term even less relevant. On the other hand, 
the range where | sin 2#k | is not small widens up with 
increasing M . Over some range of M, up to M\i = 
0.5(mj, — m)/(m y + to) < M\ on the a-pocket and up to 
M22 = 0.5(m — m x )/(m + m x ) < M 2 on the 6-pocket, 
sin2^k varies between —1 and some negative value. At 
M11 < M < Mi, sin2#k on the a-pocket varies in a 
more narrow interval, does not reach —1 and eventually 
becomes 



sin 2#k c 
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= -2 



v 77 
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= -0.98. (16) 
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The same behavior holds for sin26>k on the b— pocket 
for M22 < M < M 2 , but with 



sin26> k 



bM^hk 



= -2^ x - = -0.94. (17) 

to + m x 



III. 



GAP STRUCTURE IN THE COEXISTENCE 
STATE AS FUNCTION OF M 

From small to intermediate M: nodeless-nodal 
transition 



We now discuss in more detail the solution of the ma- 
trix equation (fi"4"f and the structure of the SC gap. We 
found in the previous section that the angular depen- 
dence of cos 20k and sin 2#k along the a and 6— pockets 
depends on M, while the solution of the matrix equation 
for gi depends on M and also on s - which, we remind, 
is the ratio of the interactions s = U1/U3. 

The solution of the matrix equation can be easily an- 
alyzed in the limits of small and large s. At small s, 33 
is small compared to 32 (33/32 = 0(s)), while 32/31 is 
given by 



32 _ N c 

31 ~ N + N c , 



(18) 



This ratio is always small because N c /N is small for 
all M (see Fig. [3]). As a result, the gaps A a (q) and A&(g) 
remain essentially constants regardless of how strong the 
SDW order is. 

In the opposite limit s» 1 the behavior is different. 
Only in a narrow range of the smallest M «l/s the gaps 
remain almost angle-independent. Outside this range 32 
is small compared to 33 (32/33 = 0(1/ s)), while 33/31 is 
given by 



53 

3i 






(19) 



Both N s and N ss become non-zero at a finite M and 
their ratio is of order one, i.e., 33 ~ g\. This leads to 
sizable angular dependencies of A a (g) and A&(g). For 
small M < Mi, N s is smaller in magnitude than N ss 
(see Fig. [3]), hence | — | < 1 and the angular variations 
of the gaps do not give rise to nodes. However, as M 
increases and approaches Mi, the ratio N s /N ss changes 
sign, becomes negative and its magnitude exceeds one 
(see Fig. [3j. The combination of this behavior and the 
fact that for M ~ Mi sin 29 reaches —1 at the tip of 
the bananas implies that the gap along the 6— pocket 
develops accidental nodes. This happens even before M 
reaches Mi because at M = Mi (when the a— pocket 
disappears), |A s /iV ss | is already larger than one, as one 
can immediately see from Eq. 1151 if one neglects in this 
equation the contributions from the a pocket. The nodes 
in Af,(q) obviously survive up to M22 because at M < 
M22 the minimal value of sin20 q remains — 1. At larger 
A/22 < M < M2 this minimal value becomes smaller 
than one, but we checked numerically that the nodes still 
survive and remain present up to M = Ma. 

To understand the gap structure in between the limits 
s <C 1 and s ^> 1, we solve the 3x3 gap equation (fT4")) 
numerically for several s with M as a running parameter, 
and for several M with s as a running parameter. In the 
two panels in Fig. [5] we plot the ratios 32/31 and 33/31 
as function of s for two representative values of M, M — 
0.09 < Mi and Mi < M = 0.12 < M 2 , and as function 
of M for s = 5. We see that \g2lg1 \ is always small, 
while 133/31! increases with increasing s and for M — 
0.12 becomes larger than 1 above a certain s. Once this 
happens, the gap on the b— pocket develops accidental 
nodes. Overall, these results show that the behavior that 
we obtained analytically at large s extends to all s > 1. 
For these values of s, the gap along the b— pocket, which 
survives up to larger M, necessarily develops accidental 
nodes around M = Ml, where the a-pocket disappears. 

In Fig. [5] we show the variations of A a (q) and A(,(q) 
for two different values of s, s = 1 and s = 5, and three 
different values of M: M < Mi, M < Mi, and Mi < 
M < M2. For s = 1, the gaps A a and A^ have no nodes 
for all M, consistent with (33/31 1 < 1 in Fig. [5] The 
only effect of the disappearance of the a— pockets is the 
switch between the maxima and the minima of the gap 
function along the remaining 6-pockets. This switch is a 
consequence of the sign change of both 32/31 and 33/31 
(see Fig. O. 

The situation is different for larger s = 5. Even for 
small M <C Mi, the gap variation along both a— and 
b— pockets become substantial. Once M becomes large 
enough to (almost) gap out the a-pockets, the nodes ap- 
pear near the tips of the 6— pocket bananas. This behav- 
ior is in full agreement with the analysis of 33/31 earlier 
in this section. 

We therefore conclude that the two conditions to ob- 
tain nodes in the superconducting gap in the coexis- 
tence phase with stripe antiferromagnetism are: (i) rel- 
atively large density-density interaction U\ (leading to 
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Figure 6: The gap structure for two values of a = U1/U3 and three different M corresponding to the cases when the a and 
6-pockets are of comparable size (left column), when the a pockets are about to disappear (middle column), and when only the 
b— pockets are left (right column). Note the appearance of gap nodes for s = 5 and M = 0.12, when one pair of reconstructed 
pockets disappear. 



magnetism) compared to the pair-hopping interaction 
C/3, and (ii) the disappearance of one of the two pairs 
of reconstructed Fermi pockets. Our analysis agrees with 
previous works that found fully gapped quasi-particle ex- 
citations in the coexistence state in the cases of small M 
(Ref. G3) and U x = (Ref. [JJ). 

The other interactions, which we listed in Sec. II but 
did not include into Hi nt in Eq. (jTTJ) , modify the value 
of s and therefore affect whether or not the nodes appear 
upon increasing M. In particular, the exchange term 



U2cU.fl 



changes s to s 



U1-U2 

Us ' 



If U2 is negative, 



it makes the development of nodes more likely, while if 
it is positive and smaller than U\, it decreases s and 
may eliminate the nodes. The inclusion of this term also 
opens up the somewhat exotic possibility of a negative s. 
Although negative s is unlikely for FeSCs, we analyzed 
the s < case for completeness and show the results in 
Fig. [Jj We see that now the behavior is non-monotonic 
with s: the nodes appear at some intermediate \s\ and 
disappear at larger \s\, already at relatively small M, 
when all four pockets are present. 



B. From intermediate to large M: nodal-nodeless 
transition 



In the previous section we showed that nodes appear 
at large enough s at intermediate values of M ~ Mi, 
where one of the two pairs of reconstructed pockets (the 
a-pocket) disappears. We now analyze what happens 
when M becomes larger than M^ and 6-pockets also dis- 
appear. 

We find that the results depend on whether we restrict 
the pairing to the FS or include into the pairing problem 
also the states which are already gapped out by SDW. 
If we restrict the pairing problem to the FS, as we did 
before, we find that the nodes in A;, are present for all 
M < A'h- At M = AI2 the remaining FS disappears 
and T c vanishes. If we do not restrict the pairing to 
the FS and take into consideration the states already 
gapped by SDW, SC persists into the region M > M 2 (see 
Refs. 14.16.24). In this last region, all states are gapped 
already above T c , and the opening of an additional SC 
gap only moves states further away from zero energy. 
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Figure 7: Gap structure for negative s and M = 0.09. As \s\ increases, the gap on the b pocket is pushed down, acquires the 
nodes, and eventually reverses sign and again becomes nodeless. For all three values of s the FS contains both a and b pockets. 



The thermal conductivity and the penetration depth in 
this last regime show exponential behavior, typical for a 
superconductor with a full gap, From this perspective, 
the evolution of the system response around M = M 2 
mimics the transformation from a nodal to a nodeless 
gap, irrespective of whether the actual SC gap by itself 
evolves around M%. 

On a more careful look, we find that the gap A& does 
change its structure near M = M 2 and becomes nodeless. 
At around the same M, T c is likely non-monotonic and 
passes through a minimum. 

The argument goes as follows. Consider first the pair- 
ing confined to the FS and solve for T c . Neglecting N c 
and N cs in Eq. (fH|) . which are small for all M and s, and 



solving for the linearized gap equation for L = log I A 
we obtain 



U: 



1 + -^ (sN ss -N)L + 



s (iV s 2 - NN SS ) L 2 = 



2 v ~~ " "'' ' V 2 

(20) 

One can verify that NN SS — N 2 is positive for all M > 

M%. For s > 1, when A& has nodes at M ~ Mi, we 

also have sN ss > N. In this situation, there is a single 

solution of (J20J) with L > 0: 



L 



1 



sU 3 (NN ss - N 2 ) 



x (y/(sN as +N) 2 -4sN2 + (sN ss - JV)) (21) 

As long as N 2 < NN SS , L w 2/U 3 N, i.e. T c does 
not depend strongly on M. However, near M — M2, 
sin26*(k;) on the b— pockets tends to a constant value c, 
see Eq. (fTT|) . (c = —0.94 for our parameters) and N 2 and 
NN SS both tend to the same value c 2 N 2 . As a result, 
L diverges as (2/U 3 s)(sN ss - N)/(NN SS - N 2 ), i.e. T c 
vanishes (see Fig [8]). One can straightforwardly show 
that L diverges at M = M 2 even if we solve the full 3x3 
linearized matrix gap equation, without neglecting the 
terms N c and N cs . 

In the same limit NNe, m N 2 we also have 



.91 



sin 26 n 






sin 29 n 



1. 



(22) 



We see that the angle-independent and the sin 29 terms 
in A5 become identical and cancel each other, i.e Af, van- 
ishes (see Eq. f|13[) ) . This is consistent with the vanishing 
oiTr. 
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Figure 8: a) Schematic behavior of T c as a function of M, if T c 
is obtained by restricting to contributions only from the FSs. 
This T c drops to zero asM-> M2. The arrow indicates the 
direction along which M increases, b) The actual T c (solid red 
line) is a superposition of contributions from the FS and from 
the gapped states (dashed lines). This T c initially decreases 
when M approaches M2 but reverses the trend and remains 
finite at M2 and even in some range of M > M 2 . The dip in T c 
is in the region where the contributions from the FS and from 
the gapped states become comparable to each other. This 
non-monotonic behavior of T c around M2 is captured by Eq. 
((24)) (see text). 



These results, however, hold only as long as we re- 
strict the integrals for N, N s , and N ss to the FS. Once 
we include the contributions from the gapped states, 
N 2 and NN SS no longer tend to the same value. The 
contributions from the gapped states contain, instead 
of L = logA/T c , the factor L D = log A/^/D 2 + T?, 
where D is the gap in the absence of SC. In a generic 
weak coupling case, these contributions would be small 
in Ljj/L compared to the contributions from the FS, but 
in our case the FS contribution to NN SS — N 2 vanishes at 
M 2 , and the contribution from the gapped states become 
the dominant ones. Similarly, contributions from the 
gapped state break the cancelation between the angle- 
independent and the sin2# q terms in A;,. 

Combining the contributions from the FS and from the 



gapped states to iV*, we have near M = M2, 



V = .V„(l+7^ 



v, = y,\c s - ls ^j- 



N„ 



N c ss + 7 S 



L 



D 



(23) 



where c s sa — c > and c ss sa c . One can introduce 
e = c ss — c 2 as a measure of deviation from M2, with 
e = at A/2. The constants 7, are all positive, at least if 
they predominantly come from gapped a— pockets. Since 
sin20 q along the a— pockets is quite close to —1 (see Fig. 
2]), we have 7 2 = 77^ to a good accuracy. Substituting 
these forms into the result for L, Eq. (f2"Tj) . and assuming 
for simplicity that s > 1, we obtain for small e 



K 2 



- K X L -K 2 = 

(7c 2 + 7;jS +27 s |c|)Z i 

JssL D > 



(24) 



where L = L(N U 3 /2), L D = L D (N Q U 3 /2). If we set 
Lo = 0, i.e., neglect the contributions from the gapped 
states, K% — 0, K\ — — c 2 , and we obtain L = c 2 /e, as 
before. If instead we set e = 0, we find that L remains 
finite and positive, and to leading order in L^ < 1 is 
given by 



Lf_—r, 



Is 



j ss +^fc 2 + 27 s |c| 



(25) 



The behavior of L at small but finite e is rather involved, 
but its main features can be understood directly from 
Eq. ([24| . Far from M2, T c is not small, implying that 
Le, is small, i.e and K\ < and K% can be neglected. 
As a result, L increases as L ~ 1/e, leading to a decrease 
in T c and, consequently, to an increase in Ljj. With in- 
creasing Ld, K\ crosses zero and changes sign. At this 
point, K 2 = 0(1), and L ~ K^/y/e ~ 1/y/t, i.e., it 
still increases with decreasing e, but more slowly then 
before. Because L €=0 ~ 0(1) given by (|25p is certainly 
smaller than l/y/e, L necessary passes through a max- 
imum near the point where K\ changes sign, and then 
decreases towards L e= $. Accordingly, T c — Ae~ 2L ' UsN ° 
passes through a minimum at some M < M2, as shown 
schematically in Fig. [8] 

Contributions from the gapped states also affect the 
interplay between the constant and the sin 2# q terms in 
At. if we include these terms, we obtain, near M = M 2 , 
instead of Eq. ([22]). 



— sin26» 
9i 






(26) 



where, we remind, in our case j s . ss > and c = —0.94. 
We see that the numerator gets smaller than one and 



the denominator gets larger than one. As a result, 
I 22 - sin20 q | becomes smaller than one, and the gap A^ 
recovers its nodeless form in the vicinity of M 2 - Very 
likely A;, remains gapless also at larger M, before T c fi- 
nally vanishes. 

We also note that our consideration is self-consistent in 
the sense that the transformation from nodal to nodeless 
gap and non- vanishing of T c at M = M2, where e = 0, 
are consistent with each other. Indeed, once the gap 
At, becomes nodeless, the pairing problem at M > M 2 
is qualitatively similar to the one for fermions with the 
dispersion 



£gapped\K ) 



D* + [v,-(k-k F )]' 



(27) 



where D is the pre-existing gap in the excitation spec- 
trum. For such model, T c is non-zero at D — (equiv- 
alent of M = M2) and extends into the region where 
D^0: 



T C (D) = T c (0)F 



D 



Tc(0) 



(28) 



where F(x) is a decreasing function of x subject to 
F(x -> 0) = 1 - x 2 and F{x) = -1/ In (1.76 - x) near 
the critical x — 1.76. 

The outcome of this analysis is that nodal SC gap ex- 
ists only in relatively narrow range of M, from M < Mi 
to M < M.2- At larger M the gap is again nodeless, 
and, furthermore, for M > M2 all low-energy states are 
gapped already above T c . The SC transition temperature 
T c is non-zero at M 2 and decreases into M > M 2 region. 
Before that, T c has a minimum roughly where the gap 
changes from nodal to nodeless. 



C. Role of the non-reconstructed pockets 

So far, we have restricted our consideration of the pair- 
ing problem to one hole and one electron pocket recon- 
structed by SDW. There are additional pockets which do 
not participate in the SDW stated. In this subsection 
we analyze to what extent these additional pockets af- 
fect our results. Specifically, we add another elliptical 
electron pocket centered at (n, 0), denote fermions near 
this pocket by dk, and solve the set of coupled linearized 
gap equations for A a (q), A(,(q), and Arf(q). The interac- 
tions involving fermions near (0, 0) and near the electron 
pocket at (jr, 0) are the same as the interactions between 
(0,0) and (0,7r), i.e., the couplings are the same Ui, U 2l 
U3, and U4. As before, we consider the model with only 
t/3 and Ui = sU 3 . We also neglect direct interaction 
between electron pockets. 

The effective pairing interactions in the a — b — d space 
are obtained by dressing the interactions by coherence 
factors associated with the transformation from c and 
f to a and b operators for the pockets at (0, 0) and at 
(0, 7r). Since the expressions for the vertices are long, we 
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refrain from presenting them. Yet, it is straightforward 
to obtain and solve the set of linearized gap equations for 
A a (q), A;,(q), and A<j(q). It is essential that the "dress- 
ing" does not involve d— fermions, hence the interaction 
has no dependence on the angle along the d— FS. As 
a consequence, A<j(q) remains angle- independent. This 
constant A^, however, affects the interplay between the 
angle-independent and the cos 2#k and sin 2#k terms for 
A a (q) and A 6 (q). 

In Fig. [HI we show the gap structure obtained from the 
solution of the pairing problem for s — 5 and the same 
three values of M as in Fig. [51 Comparing Figs. [5] and 
Fig. [9j we see that the non-reconstructed <i-pocket has 
only a minor effect on the a — b gap structure - the gap 
on the b pocket still becomes more angle-dependent with 
increasing M and develops accidental nodes at M > Mi . 
In view of this result, it is very likely that the physics 
of gap variation with M is fully captured already within 
the two-band model of one hole and one electron pocket. 



a spin at a given site towards the field and creates a 
canted two-sublattice structure, with orthogonal antifer- 
romagnetic SDW component M = M x and ferromagnetic 
component along z. 

The quadratic part of the Hamiltonian in the presence 
of the field is 



H = J2 [( £ k + <rh) <i CT c kff + (4 +Q + ah) fl +Qa f k+Q<T 
+ E M ( c kt/k + Qi + c L/k+Qt + h.c) (29) 



where h = j-IbH. We will measure h in units of 2/z, as 
with other observables. 

The transformation to the new operators is now given 
by 



aa ka 
b\«j 
where 



-a cos f/k.er sin t/k. CT 

— Sin6>k,cr — CTCOS^kcr 



Ckcr 
/k+Qcr 



(3 



IV. EFFECT OF AN EXTERNAL MAGNETIC 
FIELD 

We now consider how the SC gap structure is affected 
by an external magnetic field H. The specific goal is 
to verify whether the field increases or reduces the ten- 
dency towards the development of nodes in A;,. This 
issue is related to experiments on (Bai_ x K x )Fe2As2, par- 
ticularly to recent measurements of thermal conductivity 
in a field—. 

We direct the field along z and include into the 
Hamiltonian the Zeeman coupling between the field 
and the z— component of the spin of a fermion S z — 
(l/2)c| t CT^ /g c^. For an isotropic system, the SDW vec- 
tor M is oriented transverse to H and for definiteness 
we direct it along x. At a finite H . the system also de- 
velops a non-zero uniform magnetization, which rotates 



J 



COSflkcr 



Smtfk.er 



and 



E~ = 



e c -e f 



M 



M 2 + (Ea) 2 



ah) - 



'e c -ef 



(31) 



M 2 



Substituting the transformation into Hint in Eq. (jTTJ) 
and taking care of the fact that cos 9h,-\ and cos 9k,\. are 
now different, we obtain Eq. ([T2l but with new (7°° k , 
U™ k and U"\ given by 



U 3 



^q.k — ^q,k = — "7T (CqfC-qJ.'S'kt'S'-ki + Sq^S-^Cy^C -\^) + C^l C q -J- 5_ q ^ C k f S"- ki 
^q,k = U qM = ~~2 (CqT C '-qi C 'ktC-ki + Sqt S- q| Skf #- ki ) - ^CtfS -qlC^S-^ 



r 



where Ck, CT = cos6*k,o- and Sk,<7 = sin #k, a- The gap 
structure consistent with these IF, ,, is of the form 

q,K 

A a (q) = giCq-tC-qi+g 2 Sq- t S-qi 

+93Cqf5-qi + 9iSqfC-ql 

Ab(q) = 9lCq-fC-ql+g 2 Sq-fS-ql 

_ ff3Cqt'5-qi _ 9iSq^C-ql 

We constructed the set of coupled linearized gap equa- 



tions for gi by standard means and solved them for var- 
ious M, s = U1/U3 and h. We found that the field en- 
hances the angle variation of the gaps A a and A;,, such 
that nodes appear at smaller M and smaller s. To illus- 
trate this, in Fig. [TO] we compare the gap structure for a 
given M = 0.04 and s = 1 (same as in the upper left panel 
in Fig. [5]) without a field and with a field. We clearly see 
that the gap variations along the reconstructed FSs grow 
with the field, and for large enough field nodes appear 



11 



1.5 



I 0.5 ry> 



en 
a. 
a 
CD 



-0.5 



-a-pocket 

-b-pocket 

d-pocket 



M=0.04 



100 200 300 
length along pocket 



400 



1.5 



CD 



^v_y 



o 0.5 



-0.5 



-a-pocket 

-b-pocket 

d-pocket 



M=0.09 



100 200 300 
length along pocket 



400 



1.5 



o 0.5 



CD 



-0.5 



— b-pocket 
d-pocket 



M=0.12 



^T\ 



100 200 300 
length along pocket 



400 



Figure 9: The gap structure for s — 5 and the same values of M as used in Fig. [6] for the case when we include an additional, 
spectator pocket at (ir, 0) which does not participate in the FS reconstruction. A comparison with Fig. [5] shows that the gap 
structure is almost unaffected by the presence of the spectator pocket. 



well before M reaches the value at which the a— pockets 
disappear. The implication is that the range where the 
gap has nodes widens up with the application of a field 
and, in particular, extends to smaller M (larger dopings), 
where without a field the system was a superconductor 
with a nodeless gap. 
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Figure 10: The gap structure for h = (left panel) and h = 
0.08 (right panel) for M = 0.04 and s = 1. The application 
of the field increases angle variations of the gaps and leads 
to the appearance of nodes at a smaller M (and a smaller s) 
than without a field. 

We emphasize that the extension of the nodal region 
to smaller M is a different effect than the field-induced 
changes in M\ and M 2 . These last changes are due to 
the modification of the quasiparticle dispersion, which in 
the presence of the field becomes 4-band dispersions and 
takes the form 



K" 




(33) 



The two a-pockets are split by ±/i into a + and a~ 
pockets, and the two b— pockets are split by ±h into b + 
and b~ pockets. 

The presence of ±h under the square root splits Mi 



,m: 



and M+, 



M~ 



where M^ 2 



and M 2 into M^,iv± l 
Mi, 2 (l + 2h) and M x " 2 = Mi )2 (l - 2h), with Mi and 
M 2 given by Eqs |9] and [10] (we recall that h is mea- 
sured in units of 2fi). For h = 0.08 used in Fig. [TU] 
we get M+ = 0.118, M^ = 0.086 instead of M x = 0.102, 
and M+ = 0.205, M 2 " = 0.149 instead of M 2 = 0.177. 
Clearly, the field-induced changes in the form of recon- 
structed FSs are small compared to the changes of the SC 



gap structure associated with the field-induced change of 
the pairing interaction. The small effect of the field on 
the reconstructed FS is in agreement with previous works 
on the cuprates^. 



DISCUSSION AND CONCLUDING 
REMARKS 



We showed that the superconducting gap of an ,s H 
SC in the co-existence phase with an SDW order ac- 
quires additional angular dependence compared to the 
case when SDW order is not present. In particular, an 
isotropic pairing interaction becomes angle-dependent in- 
side the SDW state. At small SDW order parameter M, 
or when the density-density interaction between the hole 
and electron pockets is small compared to pair hopping 
interaction term, this extra dependence is weak, in agree- 
ment with previous studie o 13 ' 15 . However, when density- 
density and pair hopping interactions are comparable and 
M is large enough to gap out one out of two pairs of 
reconstructed FS pockets, the angular dependence gets 
quite strong and the SC gap on the remaining FS banana- 
like pocket develops accidental nodes near the tips of the 
bananas. At even larger M, the system bounces back 
into the fully gapped SC state which extends into the 
regime where large enough M gaps out the remaining 
pair of FS pockets. The SC transition temperature T c 
has a dip near the point where nodal SC state becomes 
nodeless. In the presence of a magnetic field, the width 
of the nodal region expands, and, in particular, the nodes 
appear at smaller M (larger doping). 

Our results offer a consistent explanation for 
the recent experimental observations^ in underdoped 
(Bai_ 2 ,K a; )Fe 2 As 2 . Performing in-plane thermal conduc- 
tivity measurements, Reid et al. found a small range of 
doping in the SC-SDW coexistence region where the ra- 
tio k/T is finite at T — 0, i.e., k is linear in T at small T. 
This linearity is generally viewed as a strong indication 
for the presence of the nodes in the SC gap. At smaller 
and larger dopings Reid et al. found that k/T vanishes 
at T = 0, as it is expected for a fully gapped supercon- 
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ductor. They also observed that the doping range where 
k/T is finite at T — expands upon the application of 
an external magnetic field. 

In terms of our model, the nodeless-nodal-nodeless 
transition observed by Reid et al. can be attributed to 
the following sequence of events: at the optimal dop- 
ing, which roughly coincides with the onset of the co- 
existence range, the SC gaps in (Bai_ x K x )Fe2As2 are 
almost isotropic, and the excitation spectrum is fully 
gapped. Decreasing x and moving to the underdoped re- 
gion in which SC and SDW start to coexist is equivalent 
to making M nonzero in our model. In this situation, 
the gaps become anisotropic, but remain nodeless over 
some range of M. Moving deeper inside the coexistence 
region is equivalent to making M larger. Eventually, M 
becomes large enough such that one of the reconstructed 
pairs of pockets vanishes. In this situation, we found (for 
s > 1) that the gap on the remaining FS pockets develops 
nodes, what leads to a finite k/T at T = 0. As the dop- 
ing level decreases even further, M continues to increase 
and the gap structure bounces back to nodeless, since 
for such pairing state T c remains non-zero even when the 
remaining pair of reconstructed pockets disappears. The 
quasiparticle spectrum becomes fully gapped again, and 
k/T vanishes at T = 0. 

In this explanation the nodeless-nodal transition is 
roughly associated with the vanishing of one pair 
of reconstructed FS pockets. An SDW-driven elec- 
tronic transition has been observed by ARPES as 
well as by transport measurements in electron-doped 
Ba(Fei_ x Co x )2As2^. In that case, however, the re- 
constructed pocket that disappears is a much smaller 
hole-like pocket associated with additional details of the 
band dispersion not captured by our simplified two-band 
model. It would be interesting to verify experimentally 
whether the transition from four to two pockets takes 
place in the hole-doped (Bai_ :c K 2 ;)Fe2As2 materials, and 
what is their relationship to the SC gap structure. Such 
transitions should have distinct signatures not only in 



the band dispersions, but also in transport coefficients. 
Furthermore, as discussed in Refs. |27||28| . the onset of 
accidental nodes should also affect the low-temperature 
behavior of several thermodynamic quantities, giving rise 
to peculiar scaling relations. Another result of our anal- 
ysis is the observation of a dip in the doping dependence 
of T c near the doping where the system undergoes the 
nodal-nodeless transition. 

Our model indeed does not include all aspects of the 
physics of the co-existence state. In particular, doping 
with holes or electrons changes the chemical potential, 
which was kept constant in our calculation. The feedback 
from this change affects the values Mi and M% at which 
reconstructed pockets disappear leaving a possibility that 
reconstructed FSs can be present even for zero doping. 
The angular dependence of the magnetic interaction may 
also preserve FS pockets even for large M — . Still, we 
believe that the physics described by our model is quite 
genetic and should hold for more realistic models. 

Finally, we point out that the disappearance of two 
out of four pockets is not the only mechanism that can 
give rise to nodes in the coexistence state. When the 
exchange interaction between unreconstructed hole and 
electron pockets is strong enough, it can induce nodes 
even in the absence of any dramatic change of the re- 
constructed FSs. Conversely, there is also a region in 
parameter space (small s) in which the coexistence with 
SDW does not generates nodes in the SC gap. 
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